CONSTRUCTING NON-COMPUTABLE JULIA SETS 

M. BRAVERMAN AND M. YAMPOLSKY 
\j^ ' Abstract. We completely characterize the conformal radii of Siegel disks in the family 

CN I corresponding to computable parameters 9. As a consequence, we constructively produce 

Vh ■ quadratic polynomials with non-computable Julia sets. 

<:^ 

t^^ i The purpose of this note is to completely characterize the conformal radii of Siegel disks 

in the family 

Pg(z) = e^^^'z + z\ 

corresponding to computable parameters 6. As one consequence, we derive the following 
statement: 
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Theorem. There exist computable complex param,eters c, such that the Julia set J^ of the 
quadratic polynomial fc{z) = z'^ + c is non- computable. 

In the end of the note we discuss the computational complexity of constructing such values 
of c, and show that it is poly-time, assuming a certain widely believed conjecture holds. 



We recall the relevant definitions of constructive analysis very briefly. The reader is referred 



1. COMPUTABILITY OF SUBSETS OF 

m 

Q \ to our earlier work [BY] for a more detailed exposition. 

^ ■ Denote by D the set of the dyadic rationals, that is, rationals of the form ^. The classical 

definition of a computable real number may be formulated as follows: 



Definition 1.1. A number x G M is computable, if there exists a Turing Machine (further 
c^ ■ abbreviated as TM) M{n) with a positive integer input n, such that for each n, this TM 

terminates and outputs a number (i„ G D with the property |a; — (i„| < 2~". 



The definition is generalized to points in M in an obvious way. 



;^ ■ Sometimes, it is desirable to use a real number as a parameter of a computation, without 

regard to its computable properties. This is achieved with the use of oracles. 

Definition 1.2. We say that : N — * D is an oracle for a real number x, if \x — 4>{n)\ < 2"" 
for all n G N. We say that a TM M*^ is an oracle machine, if at every step of the computation 
M is allowed to query the value (f){n) for any n. 

Now, for instance, we may define computable functions of a real variable (compare with the 
discussion in [Brv]). 

Definition 1.3. We say that a function / : [a,b] -^ [c,d] is computable, if there exists an 
oracle TM M'^{rn) such that if is an oracle for x G [a, 6], then on input m, M'^ outputs a 
y G © such that \y - j{x)\ < 2"™. 
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Recall that the Hausdorff metric is a metric on compact subsets of R" defined by 

(1.1) dffiX, Y) = inf{e > | X C U,{Y) and Y C U,{X)}, 

where Ue{S) is defined as the union of the set of e-balls with centers in S. 
We introduce a class C of sets which is dense in metric dn among the compact sets and 
which has a natural correspondence to binary strings. Namely C is the set of finite unions 
of dyadic balls: 

{n 
\\ B{di,ri) I where di,ri E © 

We now define the notion of computability of subsets of M" (see [Wei]). 

Definition 1.4. We say that a compact set K (^ M.^ is computable, if there exists a TM 
M{n), n eN which outputs a set Cn E C such that distji^ (C„, K) < 2~". 

It is not difficult to see that: 

Proposition 1.1. A set K (s'R'^ is computable if and only if the distance function d^ix) = 
inf{|a; — y\ \ y E K} is a computable function of a real variable. 

In this note we are concerned with computability of quadratic Julia sets Jc = J{z^ + c). 
Define the function J : C — ;► K* {K* is the set of all compact subsets of C) by J (a) = J{fc)- 
In a complete analogy to Definition 1.3 we can define 

Definition 1.5. We say that a function k, : S —^ K* for some bounded set S is computable, 
if there exits an oracle TM M'^{m) with representing x E S such that on input m, Af^ 
outputs a C E C such that dniC, k,{x)) < 2"™". 

In the case of Julia sets: 

Definition 1.6. We say that Jc is computable if the function J : d ^-^ Jd is computable on 
the set {c}. 

We have established the existence of non-computable quadratic Julia sets in [BY], based on 
the analysis of Julia sets with Siegel disks. The next chapter summarizes the relevant tools 
of Complex Dynamics. 

2. CONFORMAL RADII OF SlEGEL DISKS AND YOCCOZ'S BRJUNO FUNCTION 

Let i? : C ^ C be a rational map of the Riemann sphere. For a periodic point zq = R^{zo) 
of period p its multiplier is the quantity A = X{zo) = DW{zq). We may speak of the 
multiplier of a periodic cycle, as it is the same for all points in the cycle by the Chain Rule. 
In the case when |A| 7^ 1, the dynamics in a sufficiently small neighborhood of the cycle 
is governed by the Intermediate Value Theorem: when < |A| < 1, the cycle is attracting 
{super- attracting if A = 0), if |A| > 1 it is repelling. Both in the attracting and repelling 
cases, the dynamics can be locally linearized: 

(2.1) ^l:{R^{z)) = \-i,{z) 

where ■?/' is a conformal mapping of a small neighborhood of zq to a disk around 0. By a 
classical result of Fatou, a rational mapping has at most finitely many non-repelling periodic 
orbits. 
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In the case when A = e^'^*^, 9 E M., the simplest to study is the parabolic case when 
6 = n/m G Q, so A is a root of unity. In this case RP is not locally linearizable; it is not hard 
to see that zq G J{R). In the complementary situation, two non- vacuous possibilities are 
considered: Cremer case, when RP is not linearizable, and Siegel case, when it is. In the latter 
case, the linearizing map ip from (2.1) conjugates the dynamics of R^ on a neighborhood 
U{zq) to the irrational rotation by angle 6 (the rotation angle) on a disk around the origin. 
The maximal such neighborhood of zq is called a Siegel disk. 

Let us discuss in more detail the occurrence of Siegel disks in the quadratic family. For 
a number 9 G [0, 1) denote [ro,ri, . . . ,r„, . . .], r^ G N U {oo} its possibly finite continued 
fraction expansion: 

(2.2) [ro, ri, . . . , r„, . . .] = 



^0 



ri + 



Such an expansion is defined uniquely if and only if ^ ^ Q. In this case, the rational 
convergents Pn/Qn = [1^0, ■ ■ ■ ,^n-i] are the closest rational approximants of 9 among the 
numbers with denominators not exceeding g„. In fact, setting A = e^'^*^, we have 

IA'* - 1| > |A^" - 1| for alio < h < g„+i, h ^ g„. 

The difference |A^" — 1| lies between 2/qn+i and l-n jqn^i, therefore the rate of growth of the 
denominators g„ describes how well 9 may be approximated with rationals. 
We recall a theorem due to Brjuno (1972): 

Theorem 2.1 ([Bru]). Let R he an analytic map with a periodic point zq G C. Suppose that 
the multiplier of zq is X = e^'^*^, and 

(2.3) B{9) = J2^-^^^^^<oc. 

Then zq is a Siegel point. 

Note that a quadratic polynomial with a fixed Siegel disk with rotation angle 9 after an 
afiine change of coordinates can be written as 

(2.4) Pe{z) =z^ + e^^''z. 

In 1987 Yoccoz [Yoc] proved the following converse to Brjuno's Theorem: 

Theorem 2.2 ([Yoc]). Suppose that for 9 G [0, 1) the polynomial Pg has a Siegel point at 
the origin. Then B{9) < 00. 

The numbers satisfying (2.3) are called Brjuno numbers; the set of all Brjuno numbers will 
be denoted B. It is a full measure set which contains all Diophantine rotation numbers. In 
particular, the rotation numbers [ro, ri, . . .] of hounded type, that is with sup r^ < 00 are in B. 
The sum of the series (2.3) is called the Brjuno function. For us a different characterization 
of B will be more useful. Inductively define 9i = 9 and 9n+i = {l/6'„}. In this way, 

"n [' n— 1) ' ra? ' n+lj . . .J • 
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We define the Yoccoz's Brjuno function as 

n=l ^"^ 

One can verify tliat 

B{d) < oo^$(^) < OO. 
Tlie value of the function $ is related to the size of the Siegel disk in the following way. 

Definition 2.1. Let {U,u) be a simply-connected subdomain of C with a marked interior 
point. Consider the unique conformal isomorphism : D i-^ f/ with 0(0) = u, and 0'(O) > 0. 
The conformal radius of {U,u) is the value of the derivative r{U,u) = 0'(O). 

Let P{0) be a quadratic polynomial with a Siegel disk Ag 3 0. The conformal radius of 
the Siegel disk Aq is r{9) = r{Ao, 0). For all other 6 G [0, oo) we set r{6) = 0, and Ag = {0}. 

By the Koebe 1/4 Theorem of classical complex analysis (see e.g. [Ahl]), the radius of the 
largest Euclidean disk around u which can be inscribed in U is at least r{U,u)/4. 
We note that one has the following direct consequence of the Caratheodory Kernel Theorem 
(see e.g. [Pom]): 

Proposition 2.3. The conformal radius of a quadratic Siegel disk varies continuously with 
respect to the Hausdorff distance on Julia sets. 

Yoccoz [Yoc] has shown that the sum 

^{e) + \ogr{6) 

is bounded below independently oi9 E B. Recently, Buff and Cheritat have greatly improved 
this result by showing that: 

Theorem 2.4 ([BC]). The function 9 i-^ ^{9) + logr(^) extends to^ as a 1-periodic con- 
tinuous function. 

In [BBYl] we obtain the following result on computability of quadratic Siegel disks: 

Theorem 2.5. The following statements are equivalent: 

(I) the Julia set J{Pe) is computable by a TM with an oracle access to 9; 
(II) the conformal radius r{9) is computable by a TM with an oracle access to 9; 
(III) the inner radius inizedAg \z\ is computable by a TM with an oracle access to 9. 

We note that when 9 is not a Brjuno number, the quantities in (II) and (III) are each equal 
to zero, and the claim is simply that J{Pe) is computable in this case. 

We will make use of the following Lemma which bounds the variation of the conformal radius 
under a perturbation of the domain. It is a direct consequence of the Koebe Theorem (see 
e.g. [RZ] for a proof). 

Lemma 2.6. Let U be a simply- connected subdomain of C containing the point in the 
interior. Let V C U be a subdomain of U . Assume that dV C B^{dU). Then 



< r{U, 0) - r{V, 0) < A^/r{U,0)^/^. 
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3. Computing noble Siegel discs 

We will make use of a computability result for noble Siegel disks. The term "noble" 
is applied in the literature to rotation numbers of the form [oq, ai, • • • , a^, 1, 1, 1, . . .]. The 
noblest of all is the golden mean 7^, = [1, 1, 1, . . .]. 



Lemma 3.1. There is a Turing Machine M, which given a finite sequence of numbers 
[ttQ, ai, . . . , Cfc] computes the conformal radius r^ for the noble number 7 = [ag, . . . , a^, 1, . . .]. 

Noble (or more generally, bounded type) Sigel quadratic Julia sets may be constructed by 
means of quasiconformal surgery (cf. [Doul]) on a Blaschke product 

fJz) = e2--W^2^^_ 
■'^^ ' 1-3^ 

This map homeomorphically maps the unit circle T onto itself with a single (cubic) critical 
point at 1. The angle r(7) can be uniquely selected in such a way that the rotation number 
of the restriction p{f'y\'f) = 7- 
For each n, the points 

{!,/,(!), /,^(l),...,/,^--i(l)} 

form the n-th dynamical partition of the unit circle. We have (cf. Theorem 3.1 of [dFdM]) 
the following: 

Theorem 3.2 (Universal real a priori bound). There exists an explicit constant B > 1 
independent of 7 and n such that the following holds. Let 7 G M \ Q and n E N. Then any 
two adjacent intervals I and J of the n-th dynamical partition of f^ are B- commensurable: 

B'^\I\ < \J\ <B\I\. 

Proposition 3.3 ([He]). For each noble 7 = [oq, . . . , Ofc, 1, . . .] the Blaschke product f^ is Ki- 
quasisymmetrically conjugate to the rotation R^ : x ^-^ x + '~f mod Z. The quasisymmetric 

constant may be taken as Ki = (2maxaj) 

Let us now consider the mapping \1' which identifies the critical orbits of /^ and Pj by 

^ : /;(!) ^ p;{c,). 

We have the following (see, for example. Theorem 3.10 of [YZ]): 

Theorem 3.4 (Douady, Ghys, Herman, Shishikura). The mapping ^ extends to a 
K -quasiconformal homeomorphism of the plane C which maps the unit disk D onto the 
Siegel disk A^. The constant K may be taken as the quasiconformal dilatation of any global 
quasiconformal extension of the Ki-qs conjugacy of Proposition 3.3. In particular, K < 2Ki. 

Elementary combinatorics implies that each interval of the n-th dynamical partition con- 
tains at least two intervals of the (n -|- 2)-nd dynamical partition. This in conjunction with 
Theorem 3.2 implies that the size of an interval of the (n -|- 2)-nd dynamical partition of f^ 
is at most r" where 
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distH(fin,<9A^) <is:r", 



We now complete the proof of Lemma 3.1. Denote Wn the connected component containing 
of the domain obtained by removing from the plane a closed disk of radius 2Kt"' around 
each point of 

Qn = {P^{cy), i = 0,..., qn+2}- 
By Theorem 3.4, 

and we have 

Wn C A^ and disti:^(9A^, dW^) < e„ = 2irr". 
Any constructive algorithm for producing the Riemann mapping of a planar region (e.g. that 
of [BB]) can be used to estimate the conformal radius r{Wn,0) with precision e„. Denote 
this estimate r„. 

Elementary estimates imply that the Julia set J{P-y) C 5(0,2). By Schwarz Lemma this 
implies r{A^, 0) < 2. By Lemma 2.6 we have 

|r(A^,0)-r„|< \r{A^,0) - r{Wn,0)\ + \r{Wn,0) - rn\ < 4^ + en -^ 0, 

n— >oo 

and the proof is complete. 

4. Main result 

Definition 4.1. (cf. [Wei]) A real number r is right- computable (or right r.e.) if there exists 
a TM M{n) which outputs a non-increasing sequence of dyadic numbers Vn > r such that 
r-n \r. 

A right computable number does not have to be computable, as the definition does not 
require the existence of a constructive bound on the difference r„ — r. We will give a specific 
example of a non-computable r in §5. 

The following result gives a precise characterization of the conformal radii of quadratic 
Siegel Julia sets Jz^+c with a computable c. 

Theorem 4.1. Let r G [0,0.1] be a real number. Then r = r{c) is the conformal radius 
of a Siegel disc of the Julia set Jz^+c for some computable number c if and only if r is 
right- computable. 

Proof. We prove the "only if direction here. The remainder of the section is dedicated to 
proving the "if direction. 

We assume that c is computable, and show that r{c) is right-computable. The case r{c) = 
is obvious, we assume that r{c) > for the remainder of the proof. Recall that periodic 
orbits are dense in the Julia set Jc- Let Hn be the union of the first n repelling periodic 
orbits. Let a be the center of the Siegel disc we are considering. Note that a is a computable 
real number. 

We can algorithmically find a strictly increasing sequence {ni} C N such that C\Ui has 
a simply-connected component Wi containing a, where 

B2-(l+i){Hni) C Ui C B2-l{Hni). 

Using any constructive algorithm for computing the conformal radius [BB] we can approx- 
imate the k-th term of the sequence 

i?fc=r(52-(.-i)(W^,),a)+5-24-V 
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By Lemma 2.6, Rk -^ r{c), and moreover, {Rk} is a non-increasing sequence. Let pk be a 
dyadic approximation of Rk that we compute so that \pk — Rk\ < 2"^^. Let 

rk=pk + S- 2~K 

Then {r^} is a computable sequence of dyadic numbers. We have 

hm Tk = hm pk = hm Rk = r{c), 

k^oo fe— +00 k^oo 

and for each /c, 

rk = Pk + 3 ■ 2-'' >Rk + 2-2-'' > Rk+i + 4 ■ 2'^^^+^^ > pk+i + 3 ■ 2-(^'+i) = r^+i. 

This shows that r{c) is right-computable. D 

We now want to prove the "if direction of Theorem 4.1. Given a computable sequence 
{vn} such that r„ \ r we claim that we can construct a c such that r = r{c). It is easy to 
do in the case when r = 0, as any parabolic parameter would do. From now on, we assume 
that r > 0. We will be using the following two lemmas. The first one is Lemma 3.1 of [BY], 
and the second one is Lemma 4.2 of [BBY2]. 

Lemma 4.2. For any initial segment I = [oq, fli, . . . , a„], write u = [oq, c^i, • • • , o„, 1, 1, 1, . . .]. 
Then for any e > 0, there is an m > and an integer N such that if we write j3 = 
[qq, Oi, . . . , a„, 1, 1, . . . , 1, A^, 1, 1, . . .], where the N is located in the n + m-th position, then 

<^{u)+e < $(/5) < $(cj) + 2e. 

Lemma 4.3. For u as above, for any e > there is an rriQ > 0, which can be computed from 
(oo, ai, . . . , a„) and e, such that for any m > mo, and for any tail I = [an+m, o,n+m+i, ■ ■ ■] if 
we denote 

jj = [Oi, 02, . . . , (In, 1) I) • • • 5 I) C^ra+nn 0,n+m+lj • • -J; 

then 

$(/30 >$H-e. 

Using Lemma 3.1, we can get a computable version of Lemmas 4.2 and 4.3. 

Lemma 4.4. For any given initial segment I = [ao, ai, . . . , a„] and rriQ > 0, write u = 
[ao, Oi, . . . , a„, 1, 1, 1, . . .]. Then for any e > 0, we can uniformly compute m > mo, an 
integer t and an integer N such that if we write /3 = [ao, ai, . . . , a„, 1, 1, . . . , 1, N, 1,1,...], 
where the N is located in the n + m-th position, we have 

(4.1) r{u)-2e <r{(3) <r{uj)-e, 

(4.2) $(/3) > $(cu), 
and for any 

7 = ['^Oj 0,1, • • • 5 0,n, I5 I5 • • • ) 1) -^) 1) • • • 5 I5 C„+m+t+lj Cn+m+t+2, • • •], 

(4.3) $(7) > $(cj) - 2-". 
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Proof. We first show that such m and A^ exist, and then give an algorithm to compute 
them. By Lemma 4.2 we can increase ^{uj) by any controlled amount by modifying one 
term arbitrarily far in the expansion. 

By Theorem 2.4, f : 6 \-^ '^(^) + logr(^) extends to a continuous function. Hence for any 
Eq there is a 5 such that \f{x) — f{y)\ < Bq whenever |a; — 1/| < S. In particular, there is an 
mi such that |/(/5) — f{uj)\ < Eq whenever m>mi. 

This means that if we choose m large enough, a controlled increase of $ closely corresponds 
to a controlled drop of r by a corresponding amount, hence there are m > ttiq and N such 
that (4.1) holds. (4.2) is satisfied almost automatically. The only problem is to computably 
find such m and N. 

To this end, we apply Lemma 3.1. It implies that for any specific m and A^ we can compute 
r{(3). This means that we can find the suitable m and A^, by enumerating all the pairs (m, A^) 
and exhaustively checking (4.1) and (4.2) for all of them. We know that eventually we will 
find a pair for which (4.1) and (4.2) hold. 

Finally, t exists and can be computed by Lemma 4.3. D 

We are now ready to give an algorithm for computing a rotation number 6 for which 
r{6) = lim \ r„. c is easily computed from 6. The algorithm works as follows. On stage k 
it produces a finite initial segment Ik = [oq, . . . , dnik] such that the following properties are 
maintained: 

(1) /o=[]; 

(2) Ik has at least k terms, i.e. rrik > k; 

(3) for each k, Ik+i is an extension of Ik] 

(4) for each k, denote 7^ = [4, 1, 1, . . .], then r^ + 2"^''+^) < 7fc < ^fc + S"*"; 

(5) for each k, $(7fc) > $(7^-1); 

(6) for each /c, for any extension 

/5 = [Ik-, bnik+ly bmk+2, • • •], 

$(/3) > $(7fc) - 2-K 

The first three properties are very easy to assure. The last three are maintained using 
Lemma 4.4. By this Lemma we can decrease r{'~fk-i) by any given amount (possibly in more 
than one step) by extending Ik^i to J^. Here we use the facts that the r^'s are computable 
and non- increasing. 

Denote 6 = limfe^oo7A:- The continued fraction expansion of 6 is the limit of the initial 
segments Ik- This algorithm gives us at least one term of the continued fraction expansion 
of 6 per iteration, hence we would need at most 0{n) iterations to compute 6 with precision 
2^" (in fact, much fewer iterations would suffice). It remains to prove that, in fact, 6 is the 
rotation number we are looking for. 

Lemma 4.5. The following equalities hold: 

$(^) = lim $(7fc) and r{6) = lim r{'~fk) = r. 

fc— >oo k-^00 

Proof. By the construction, the limit 6 = lim.'-fk exists. We also know that the sequence 
f^ilk) converges to the number r = lim \ r^, and that the sequence $(7a;) is monotone 
non-decreasing, and hence converges to a value ip {a priori we could have ip = 00). By the 
Caratheodory Kernel Theorem, we have r{6) > r > 0, so ^{6) < 00. On the other hand, by 
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the property we have maintained through the construction, we know that $(^) > $(7^) — 2"^^ 
for all k. Hence $(^) > ip- In particular, ip < 00. 
From [BC] we know that 

(4.4) V^ + logr = lim(<l>(7fc) + logr(7fc)) = $(6') + log r(6'). 

Hence we must have ip = $(^), and r = r{6), which completes the proof. D 

5. Constructing non-computable Julia sets 

Theorem 4.1 provides us with a tool for constructing explicit parameters c for which Jz-i+c 
is non-computable. In fact, for any right-computable r which is not a computable number, 
we obtain such a set. 

In particular, we can give an explicit construction of a Julia set parameter, for which 
computing the Julia set with an arbitrarily high precision would allow us to solve the halting 
problem. 

Let R{x, t) be the predicate which evaluates to 1 if and only if x is a valid encoding of 
a Turing Machine which terminates after exactly t steps. Then for any fixed x, R{x, •) 
evaluates to 1 at most once, and 

H{x) = 3t R{x,t) 

is the non-computable Halting predicate. Let {r^} be the following computable sequence of 
dyadic numbers 



^^ = 17^ i-EE4-^-^(^'^) 



x=l t=l 

Tj: is obviously a non- increasing sequence, and it converges to 

1 / °° 
r = lim rfe = — 1 - ^^ 4~^ ■ H{x) 

k^oo 16 \ -^^ 

\ x=l 

Hence r is a right-computable number, and there is a computable parameter c such that 
r{c) = r. Now note that computing r = r{c) with precision 4"'^^'^^) would require evaluating 
H[l), H{2), . . . , H{x), which is impossible. 

6. The complexity of ^'s for which J^ is non-computable 

In this note we have demonstrated that one can explicitly compute a parameter 6 such that 
the corresponding Julia set Je = Jp^ is non-computable. In fact, computing Jg would allow 
us to solve the Halting Problem. It is a natural question to ask whether such a parameter 6 
can be computed fast. That is, how much time is required to produce a 2~"-approximation of 
6. The algorithm for generating 6 presented above does not yield any explicit bound on the 
computational complexity of 6. In this section we show that under a reasonable conjecture 
6 can be shown to be poly-time computable. That is, there is a polynomial p{n) such that it 
takes at most p{n) steps to produce a 2~^ approximation of 9. The conjecture is a stronger 
version of Theorem 2.4. 

Conjecture 6.1. The 1-periodic continuous function f : 6 h^ ^{6) + logr{6) has a com- 
putable modulus of continuity. In other words, there is a computable function yU : N ^ N 
such that 1/(^1) - /(^2)| < 2-" whenever \9i - ^2! < 2-^("). 
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Note that in Conjecture 6.1 we do not put any restrictions on the growth of fi besides it 
being computable. fi{n) can grow extremely fast, as some computable integer functions do. 
In fact, the following much stronger conjecture has been put forward by Marmi, Moussa, 
and Yoccoz in [MMY], and is generally supported by the available evidence: 

Marmi-Moussa-Yoccoz Conjecture. [MMY] The function f : 9 \-^ $(^) + logr(6') is 
Holder of exponent 1/2. 

Assuming the weaker Conjecture 6.1, we show that the ^'s from Theorem 4.1 (and hence the 
c's) can be made poly-time computable. Of course, the strengthening of Theorem 4.1 only 
comes in the "if direction. We only state that direction here. 

Theorem 6.2 (Conditional). Suppose Conjecture 6.1 holds, and suppose there is a com- 
putable sequence ri,r2, ■ ■ ■ of dyadic numbers such that 

• {vi} is non-increasing, ti > r2 > . . ., and 

• limj^oo^ = r. 

Then there is a poly-time computable 9 (and hence a poly-time computable c = c{9)) such 
that r{c) = r. 

Proof. The proof goes along the lines of the proof of the "if direction of Theorem 4.1. We 
outline the modifications made to the proof here and leave the details to the reader. The 
key difference is that in the proof of Theorem 4.1 we used Lemma 4.4 to perform a step 
in decreasing the conformal radius from r{'~fk-i) to r{'~fk)- The algorithm there is basically 
an exhaustive search, which of course could take much more than a polynomial time in 
the precision of 7^ to compute. By assuming Conjecture 6.1 we can deal with <l>(7fc_i) and 
$(7fe) instead of the r(«)'s. We have an explicit formula for $ that converges well, and we 
can compute the continued fractions coefficient to make $(7^) close to whatever we want 
relatively fast. 

The step of going from 7^-1 to 7^ is as follows. First, we do the following computations: 

• compute dk which is the "drop" in r we are trying to achieve; we want dk/2 < 
log(r(7fc_i)) - log(r(7fc)) < 4; 

• compute using the function /x a value 6k such that \f{x) — f{y)\ < dk/S whenever 
\x -y\ < 6k. 

We have no a priori bound on how long these computations would take, but we would still 
like to be computing 9 in polynomial time. To achieve this, we use I's in the continued 
fraction expansion of 9 to "pad" the computation. 

When asked about the value of 9 with precision 2^" which is higher than what the known 
terms of the expansion [Ik-i] can provide, we do the following: 

• try to compute dk and 6k as above, but run the computation for at most n steps; 

• if the computation does not terminate, output an answer consistent with the initial 
segment [4_i, 1, 1, . . .,1 ]; 

• if the computation terminates in less than n steps proceed as described below. 

Note that so far the computation is polynomial in n. For some sufficiently large n the 
computation will terminate in n steps, at which point we will have computed dk and 6k. We 
then add more I's to the initial segment to assure that \'~fk-i — lk\ < 6k. 
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Recall that our goal is to assure that 

4/2 < log(r(7fc_i)) - log(r(7fc)) < 4. 

Withe the current initial segment for 7^ we have \'~fk-i ^Ikl < ^k, and hence in the difference 

log(r(7fe_i)) - log(r(7fe)) = $(7^) - $(7^-1) + (/(7fe-i) - filk)) 

the last term is bounded by 4/8. This means that for the current step it suffices to increase 
<l>(7fe) relative to $(7fe_i) by between |4 and |4- 

Let M be the total length of Ik-i and the I's we have added, and let us extend the 
continued fraction by putting iV G N in the M + 1-st term, and all I's further. Increasing 
M if necessary, we can ensure an approximate equality 

<^{-fk) ^ Hlk-i) + a{N) log N 

up to an error of ^c?fc. Let Pm/Qm be the M-th convergent of the resulting continued fraction. 
Recall that on an input n we need to compute 6 with precision 2"" in time polynomial in n. 
If 2^" > I/^qm, then we do not need to know anything about A^ to compute the required 
approximation. Suppose 2~" < I/^qm, which means n > loggM/2. And we have time 
polynomial in logg^ to perform the computation. 

Note that M = 0{logqM)- It is also not hard to see that a{N) < 2^/^, so in order to 
have a change by ~ 3(^^/4 we must have A^ > e^*-^ -*, hence by making M sufficiently large 
(depending on the value of rffc), we can guarantee that N > e^ . This means that we can 
approximate a{N) with the truncated function $ at the M-th convergent of the continued 
fraction. Write Pm/qm = [(11,0,2, ■■■ , om], and denote 

P = [Oi, 02, ... , Om] ■ [02, ^3; • • • ; ^m] ■ ■ ■ ■ ■ [OM-I) ^m] ' [Om]- 

Then (3 approximates a{N) within a very small relative error. In particular, we can assure 
that 



32J \ 32^ 

In time polynomial in logg^ we can compute the exact expression for j3 using rational 
arithmetic: f3 = p/q. Now we can estimate A^ and write it as e^'^''^^^ in time polynomial in 
log(gM). From there we can continue by adding enough I's to get Ik and 7^ = [Ik, 1,1,.. .]. 
By the construction, it would give us the necessary decrease in the value of ^(7^). D 
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